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ABSTRACT 

Aims. Our goal is to investigate the X-ray timing and spectral variability of the high-mass X-ray binary 1A 1 118-615 during a type-II 
outburst. 

Methods. We performed a detailed color, spectral and timing analysis of a giant outburst from 1A 1 118-615, using RXTE data. 
Results. We report the discovery of a variable quasi-periodic oscillation (QPO) in the power spectral density of 1A 1 118-615, with a 
centroid frequency of ~0.08 Hz. The centroid frequency of the QPO correlates with the X-ray flux, as expected according to the most 
accredited models for QPO production. For energies above ~4 keV, the QPO rms variability decreases as the energy increases. Pulse 
profiles display energy dependence, with a two-peak profile at lower energies, and a single peak at higher energies. From spectral 
analysis, we confirm the presence of a cyclotron absorption feature at ~60 keV, the highest value measured for an X-ray pulsar. We 
find that the spectral parameters (photon index, cutoff energy, iron fluorescence line strength) display a marked dependence with 
flux. We detect two different levels of neutral hydrogen column density, possibly due to the Be companion activity. We report 
for the first time a correlation between the timing and spectral parameters in an X-ray pulsar. All the correlations found between 
spectral/timing parameters and X-ray flux are present up to a flux of ~6xl0~ 9 erg cm -2 s -1 , when a saturation level is reached. We 
propose that the saturation observed corresponds to the minimum extent of the neutron star magnetosphere. We estimate the magnetic 
field of the neutron star from two independent ways, using results from spectral (cyclotron line energy) and timing (QPO frequency) 
analysis, obtaining consistent values, of ~7 - 8 x 10 12 G. Results from the comprehensive spectral and timing analysis are discussed 
in comparison with other X-ray pulsars. 
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1. Introduction 

Be/X-ray binary systems (Be/XRBs) constitute a subclass of 
neutron-star high-mass X-ray binaries (HMXBs) in which a 
neutron star orbits a Be star. In these systems, the Be-stars are 
deep inside their Roche lobes, as is indicated by their generally 
long orbital periods and by the absence of X-ray eclipses and of 
ellipsoidal light variations. The X-ray emission from Be/XRBs 
tends to be extremely variable, ranging from quiescence, to 
periodic, short (Type-I) outbursts separated by the orbital period 
(Lx ~ 10 36 - 10 37 erg s _1 ), to giant (Type-II) transient outbursts 
lasting weeks to months (Lx > 10 37 erg s _1 ). Be/XRBs show 
hard X-ray spectra which, in combination with the regular X-ray 
pulsations, indicate that the compact object m ust be a strongly 
magnetized neutron star (see Negueruela 20O7I for a review). 

1A 1118-615 is an accretion-powered X-ray pulsar which 
was discovered i n December 1974 by the Ariel V satellite 
(Evles et a O 19751) during an observation of Cen X-3. The same 
data r evealed pulsations with a period of 405 ±0.6 s dives et al.l 
fl97l . The X-ray spectrum was fitted by a highly-absorbed 
(Nh ~ 6 x 10 22 crrT 2 ) power law with photon index ~ 1. The 
remarkable measured absorption indicated that the circumstellar 
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environm ent played a s ignificant role in absorbing the X-ray 
emission (lives et al.lll975l) . In 1979 and 1985, the system was 
observed by Einst ein and EXOSAT satellites, respectively 
dMotch et alJ Il988l) . Only weak signal was detected in both 
occasions, showing that low-level accretion was occurring. No 
spectral or temporal analysis could be performed. 

The second giant outburst from 1A 1118-615 since its 
discovery was detected in January 1992 after a long period 
of quiesc e nce. A multiwavelength study was performed by 
ICoe et all d!994l) . They included BATSE data in the en- 
ergy range 20-100 keV, supported by optical, IR and UV 
observations. The pulsed X-ray spectrum was fitted by a single- 
temperature, optically-thin, thermal bremsstrahlung model 
with a temperature of 15 keV. A spin-up of 0.016 s day -1 was 
observed during the decay of the outburst. Lower frequency 
observations showed very strong Ha emission (EWh„= -80A), 
and high IR excess from the disk of the companion star 
(E(J - A") =0.47), indicating a very large circumstellar disk 
dCoeet ai.lil99"l . 

In cont rast to what had been previously found by 

iMotch et all d 19881). pulsations from 1A 1118-615 were 
detected bv lRutledge et all d2007l) at low luminosity (L x S 10 35 
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erg s 1 ), making it the third HMXB to exhibit pulsed emission 
at low luminosities , after A 0535+26 dNegueruela et alJ [2000) 
and 4U 1145-619 dMereghetti et al.lll987l) . After 17 years of 
quiescence, a third major outburst was observed at the beginning 
of 2009. Using data from t he Ro ssi X-ray Timing Explorer 
(RXTE), DoroshenkoetalJ d201Q|) measured the broad-band 
(3-5-120 keV) spectrum of 1A 1 118-615 for the first time. From 
observations corresponding to the peak of the outburst, they 
detected a cyclotron resonant scattering feature (CRSF) at ~55 
keV and the presence of an iron emission line at 6.4 keV. 

1A 1118-615 is associated with the Q 9 5 Ill-Ve star He3- 
640, w i th an estimated dista nce of 5+2 kpc dJanot-Pacheco et al] 
119811) . IVillada et all dl999h detected strong spectral variability 
of He3-640, with a correlation between the hard X-ray flux 
and the Ha equivalent width. They proposed the presence of 
extended, expanding material around the Be star that sometimes 
becomes dense enough in the vicinity of the neutron star to 
form an accretion disk; this increases the mass accretion rate, 
hence the X-ray flux. The formation of an accretion disk during 
Type II outbursts is still under debate. However, t he discovery 
of qu asi-periodic oscillations in some systems (see iJames et alJ 
2010, and references therein) would support this scenario. The 
presence of an accretion disk also helps explain the large and 
steady spin-up rates seen during the giant outbursts, which 
are diffic ult to account f or by means of direct accretion. More 
recently, iLin et al.1 (12 010) performed a Swift spectral analysis of 
1A 1118-615 during the 2009 outburst and a flaring event two 
months later and found evidence for emission from an accretion 
disk during the flare event, suggesting that the polar cap and the 
accretion disk emission might co-exist after an outburst. 

The orbital period of the binary is still unknown, but it 
can be estimated in the range o f 400-800 d ays from the pulse 
period/orbital period diagram of Corbet] d!986l) . 

In this paper, we present a detailed spectral and temporal 
study of the third major outburst of the system using RXTE data. 
Unlike previous work, we focus on the aperiodic variability of 
the source, reporting the discovery of a ~0.08 Hz QPO, we study 
the evolution of the spectral and timing parameters through- 
out the entire outburst, and perform pulse-phase spectroscopy. 
Section 2 describes the observations and gives information on 
the data reduction. In Sect. 3 we report the results obtained. In 
Sect. 4 we interpret our findings and compare the variability ob- 
served inlA1118-615 with that reported in other HMXBs. 

2. Observations and data analysis 

The increasing flux from 1A 1118-615 at the beginning of 
January 2009 triggered a public RXTE Target of Opportunity 
(ToO) campaign of observations. TableQ]shows the observation 
log. 



Table 1. Journal of RXTE observations. 



N. of 


Proposal 


MJD 


On-source 


pointings 


ID 


range 


time (ks) 


1 


94412 


54838.4 


1.30 


25 


94032 


54841.2-54865.0 


8.90 



In this work, data from both the Proportional Counter Array 
(PCA) and the High Energy Tim ing Experiment (HEXTE) 
onboard RXTE dBradt et all 1 19931) were employed. The PCA 
consists of five Proportional Counter Units (PCUs) with a total 
collecting area of ~6250 cm 2 , operates in the 2-60 keV range 
and has a nominal energy resolution of 18% at 6 keV. The 
HEXTE is constituted by 2 clusters of 4 Nal/Csl scintillation 
counters, with a total collecting area of 2 x 800 cm 2 , sensitive in 
the 15-250 keV band with a nominal energy resolution of 15% 
at 60 keV. Both instruments have a maximum time resolution of 

Data reduction was performed using HEASOFT version 
6.9. An energy spectrum was obtained for each one of the 26 
observations, after filtering out unsuitable data according to the 
recommended criteria 1 , employing Standard 2 mode data from 
the PCA (PCU2 only) and Standard (archive) mode from the 
HEXTE. At the time of the observations, the HEXTE Cluster A 
had ceased its modulation between on-source and off-source po- 
sitions, while Cluster B was still rocking: for this reason we only 
employed Cluster B data, excluding detector 2 since its inability 
to measure spectral information. The PCA and HEXTE spectra 
were extracted, background subtracted and deadtime corrected. 
For the PCA, the 3-30 keV energy range was retained, while 
the HEXTE provided a partially overlapping extension from 25 
to 100 keV. For each observation, the two resulti ng spectra wer e 
simultaneously fitted with XSPEC v. 12.6.0 dArnaudl fl996). 
During the fitting, a systematic error of 0.6% was added to the 
spectra. 

Power spectral density (PSD) was computed using PCA 
Good_Xenon data. In one case (observation 94412-01-01-00, 
the first of the data set) Event mode data were employed. We 
first extracted, for each observation, a light curve in the en- 
ergy range 2-20 keV (channels 0-49) with a time resolution 
of 2~ 5 s. The light curve was then divided into 128-s segments 
and a Fast Fourier Transform was computed for each segment. 
The final PSD was computed as the average of all the power 
spectra obtained for each segment. These averaged power spec- 
tra were logarithmically rebinned in frequency and corrected 
for dead time effec ts according to the prescriptions given in 
iNowak et alJ (1 19991) . Power spectra were normalized such that 
the integral over the PSD is equal to the squared fractional 
rms amplitude, according to the so-called rms -normalization 
dBelloni & Hasingerfll990tlMivamoto etafll 19911) . 

3. Results 

1A 1118-615 went into outburst at the beginning of 2009. A 
flux clearly above the background was detected on MJD 54840. 
The total duration of the outburst was 27 days. A maximum 
3-30 keV X-ray luminosity of 2.8 x 10 37 erg s _1 was measured 
at MJD ~ 54845.3, assuming a distance of 5 kpc. 

Figure[T]shows the PCA light curve and color behavior dur- 
ing the outburst, as directly derived from PCU2 count-rate. Each 
point in the diagrams corresponds to a RXTE pointing. The color 
definition is the following: soft color (SC): 7-10 keV / 4-7 
keV; hard color (HC): 15-30 keV / 10-15 keV. Different sym- 
bols mark the different phases of the outburst, the rising (open 

1 Among which, elevation ftrom the Earth greater than 
10° and pointing offset lower than 0.02°; see PCA digest at 
http://heasarc.gsfc.nasa.gov/docs/xte/pca_news.html 
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Fig. 1. PCA light curve and color behavior during the 2009 outburst. 
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Fig. 2. Example of low-luminosity (Lx = 2.3 1 x 10 36 erg s , left panel) 
and high-luminosity PSD (L x = 2.50 x 10 37 erg s -1 , right panel) with 
fitted components: a power law (dotted line), L\ (solid line), Li (dash- 
dotted line) and the QPO (dashed line). On the top of each panel the 
RXTE observation ID is reported. 



peak. 

In order to have a unified phenomenological description of 
the timing features during the outburst, we fitted all the noise 
components, with the exception of the poorly constrained lower- 
frequency one, with Lorentzian functions, each one denoted as 
Li, where i determines the type of component. The characteristic 
frequency v max of L, is denoted v,. This is the frequency where 
the component contributes most of its variance per l ogarithmic 
frequency interval and is defined (iBelloni et a l. 2002) as 



yjv 2 Q + (FWHM/2) 2 



(1) 



where vo is the centroid frequency and FWHM the full width at 
half maximum of the Lorentzian component. All the frequencies 
reported in this work refer to characteristic frequencies v max . To 
fit the noise below ~0.05 Hz, a power law was employed. 



3.1.1. Discovery of a QPO 

A QPO is detected in the PSD of 1A 1118-615 during the 
brightest observations. When the X-ray luminosity goes below 
~40% of the maximum luminosity, that is for flux below 
4xl0~ 9 ergs cirT 2 s _1 , the QPO disappears. Figure|2]shows two 
examples of PSD at different luminosity to allow comparison: 
the QPO is evident in the high-luminosity plot (right panel of 
Fig. |2), while it is undetectable in the low-luminosity one (left 
panel of Fig. |2). In the figure, the broad-noise components are 
also shown. 

The QPO was fitted by a narrow Lorentzian. Its maximum 
frequency, vqpo, is ~0.08 Hz, and the measured quality factor 
Q, obtained as Q — vq/FWHM, varies in the range 3-9. Given 
the relatively large Q value, the centroid frequency vo coincides 
approximately with the maximum frequency vqpo Figure [3] 
shows the trend followed by the characteristic frequency and the 
rms of the QPO feature during the outburst. Both parameters 
fall within the typical range seen in HMXBs. The flux shown 
in Fig. [3] is calculated in the 3-30 keV range from the best-fit 
spectral model (see Sect. 13. 2t 

We have also studied the energy dependence of the QPO 
feature. We extracted PSD in the following energy bands, 2-4, 
4-7, 7-10, 10-15, and 15-30 keV, using data from the PCA. 
We found evidence for a decrease of the QPO variability with 
energy, for energies higher than 4 keV. Similar results are ob- 
tained repeating the analysis with data from different observa- 
tions. Results are shown in Fig.|4]for observation ID 94032-04- 
02-02. 



squares), the peak (filled circles) and the decay (filled triangles). 
When not shown, the errors are of the same size as the points. 
As one can see, both colors correlate with flux. 

3.1. Timing analysis 

We obtained one average PSD for each observational interval 
(pointing) as explained above. Our selection of time resolution 
and length of segments implies a frequency range of ~0. 008-16 
Hz. Thus, peaks derived from the neutron star's pulsations do 
not affect the analysis of the PSD due to their very low fre- 
quency, being v pu i se ~ 2 mHz the frequency of the fundamental 



3.1.2. Broad-noise components 

The power spectra continuum was fitted with the sum of, at 
most, three components. The low-frequency noise (<0.05 Hz) 
is accounted for by a power law, which is not well constrained. 
Above 0.05 Hz, the noise is well described by two zero-centered 
Lorentzians. L\ dominates below 1 Hz, with a characteristic 
frequency varying between 0.03-0.28 Hz. When L2 is not 
present, it also describes the noise up to 16 Hz. This component, 
together with the power law, is the only one that is constantly 
present along the outburst. With the highest fractional rms (up 
to 28%), this is the best constrained component. L2 accounts 
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Fig. 3. QPO parameters vs. calculated 3-30 keV flux: characteristic fre- 
quency (upper panel) and fractional rms (lower panel). Different sym- 
bols mark the different phases of the outburst: rising (empty squares), 
peak (full circles), and decay (full triangles). 



K 6 



C 4 







< 


» ' 



10 20 
energy (keV) 



30 



Fig. 4. Fractional rms variability for the 0.08 Hz QPO as a function of 
energy. 



instead for the high-frequency noise, with V2 ~ 0.5-2.1 Hz. 
It is not detectable in the last four pointings due to the low 
signal-to-noise in the corresponding frequency range at low 
luminosities. Its rms varies between 4%-14%. 

In Fig. |5] we present the behavior of the broad-noise com- 
ponents L[ and L2 as a function of the measured flux, while in 
Fig.|6]the relation between the characteristic frequencies v\ and 
V2 is shown. As the flux increases, the characteristic frequency 
increases and the amplitude of variability decreases. This trend 
is significant for L\, while is much weaker for Li. Another inter- 
esting feature in Fig.[5]is that above a certain flux the variation 
of the timing parameters with X-ray flux saturates. Due to the 
observational gap it is difficult to pin down the exact value of 
the flux at which this flattening occurs, but it lies in the range 
6.5-7.5xl(T 9 ergcirr 2 s _1 . 

3.1.3. X-ray pulsations 

The 407-s X-ray pulsation is clearly seen in the light curves 
of 1A 1118-615 at all energies. However, the precise de- 
termination of the spin period is hampered by its relatively 
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Fig. 5. Frequency and fractional amplitude variability as a function of 
X-ray flux for the L, (first two panels from above) and L 2 components 
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ent phases of the outburst: rising (empty squares), peak (full circles), 
and decay (full triangles). 
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Fig. 6. Relation between the characteristic frequency of the two broad- 
noise components L\ and Li. Different symbols mark the different 
phases of the outburst: rising (empty squares), peak (full circles), and 
decay (full triangles). 



long value (~407 s) in comparison with the duration of the 
observational intervals. We selected one of the observations 
with the longest on-source time (MJD 54843.9750) and derived 
the spin period by epoch-folding, that is, by folding data over 
a period range and obtaining the chi-square of the folded light 
curve. The period that gives the highest x 2 1S taken as the 
best-fit period. We found 407.5 s. This period was refined by 
cross-correlating several pulse profiles obtained by cutting the 
light curve into 1230-s segments with a template made up from 
the observation-averaged profile. The phase shifts were fitted 
to a linear function, whose slope represents the correction in 



E. Nespoli and P. Reig: Discovery of a QPO in 1 A 1 118-615: correlated spectral and aperiodic variability 



5 



ou 
40 


-_' ' 1 1 1 1 ' 


1 i i i 

2 


1 

-4 


keV ; 


100 
80 


30 


f29% 


J 


1 




60 
40 


140 

1 *Pfl 

L(C\J 


: i i i i i ti 


III 

4 


t+ 
-7 


Tl 1 ; 
keV^ 


100 
80 


100 


^ 








60 


80 


^35% 

:l 1 1 1 1 1 


III 


++ 




40 


120 




7- 


10 


keVj 


20 


100 










15 


80 


: 30% VT^ 








10 

5 


160 
140 


J 1 1 1 1 1 


III 
10- 


15 


keVJ 


8 


120 










6 


100 










4 


80 


^33% 





IV 




2 



^39% \ 


i i i | i i i i | i i i i _ 
15-30 keV^ 


J 1 1 1 | 
^45% , 


1 1 1 1 1 1 1 

30-£ 


i i i i - 
50 keV~ 

MM" 


_ 1 1 1 1 | 
z 53% 


1 1 1 1 1 1 1 
60-1 


30 ke\t 

MM" 


_ 1 1 1 1 | 

^64% 

- 


1 1 1 1 1 1 1 
100-2E 

!V, , i , , , , 


keVJ 



0.5 1 1.5 

Pulse phase 



0.5 1 1.5 
Pulse phase 



Fig. 7. Folded light curves of 1A 1118-615 in five energy bands. The 
folding epoch is MJD 54843. Profiles with energy below 30 keV were 
obtained from PCA light curves, while above this energy from HEXTE 
data. Two complete cycles are shown for clarity The number in the 
lower-left corner represents the pulse fraction. 



frequency to the actual period. A new template was derived for 
this corrected period, and the process repeated. The error was 
estimated from the uncertainty of the slope from the linear fit. 
The final spin period was 407.77+0.08 s. Pulse profiles were 
obtained by folding light curves at different energies with the 
best-fit period. 

Figure Q shows the pulse profile at various energy bands. 
We find a significant change of the pulse profile with energy, 
with the pulse profile showing a double peak structure below 10 
keV. The intensity of the second peak decreases as the energy 
increases. Above lOkeV the second peak completely disappears 
and the pulse profile becomes nearly sinusoidal. The pulse frac- 
tion, defined as (/ max - / m i n ) / (/ max + AninX increases as the energy 
increases above 10 keV. Below 10 keV there is a relative mini- 
mum at energies 4-7 keV , in good agreement with the results of 
IDoroshenko et al.l d2010l) . 

3.2. Spectral analysis 

Energy spectra were best fitted with an absorbed power law 
with exponential high-energy cutoff and a gaussian line profile 
at ~6.5 keV to account for Fe Ka fluorescence. The inclusion 
of a cyclotron absorption feature with central energy of ~60 
keV was a necessary step in order to obtain acceptable fits. The 
obtained reduced chi-square was x^ ed = 0.9-1.6, for 72 degrees 
of freedom (d.o.f.). 

Figure [8] shows the behavior followed by the main spectral 
parameters as a function of the flux. Along the outburst, the pho- 
ton index anti-correlates with flux, implying, as also the colors 
show, softer emission at low flux and harder emission at high 
flux. A similar trend is followed by the cutoff energy, with lower 
values in correspondence to higher flux. The best-fit hydrogen 
column density, although it can be only poorly constrained by 
the PCA, remains approximately constant during the outburst 
evolution at ~ 3 x 10 22 cm" 2 ; only a few observations at the 



very end of the available data sample and one observation at 
the beginning show higher values: in these cases, very poor fits 
were obtained by fixing the value of A^h at ~3xl0 22 cm" 2 , with 
X^ ed ~2-3 (75 d.o.f.). By letting this parameter free to vary, the 
retrieved ^ 2 ed lowered to ~1.05 (76 d.o.f.), and the value of A^h 
was ~7xl0 22 cm" 2 . 

Using the so urce infrared excess E(J - K) = 0.47 
(ICoe et al.l [T99l and assuming the mean extinction law 
(Rv = 3.1), we obtained the total measur ed visual extinctio n 
from A V /E(J - K) = 5.82 ± 0.10 dRieke & Lebofskvl [T985). 
We converted it to the corresponding hydrogen column density 
value from N u /A v = 1 .79 ± 0.03 x 10 21 atoms cm 2 mag 
dPredehl & Schmittl 1995) in order to compare the interstellar 
value of A^h, as obtained from the extinction, with the total 
one, provided by X-ray data. We obtained a visual extinction of 
Ay = 2.7 mag, which corresponds to A^h = 4.9 x 10 21 atoms 
cm -2 . This value is one order of magnitude lower than the one 
obtained by our X-ray data. This difference can be explained 
assuming that the source of extinction measured by the X-ray 
band differs from that measured in the optical and infrared 
bands. The derived optical/IR A^h measures the interstellar 
extinction, while the A^h obtained from the X-ray spectrum is 
also affected by absorbing material locally, that is, in the vicinity 
of the compact object. An additional local source of extinction 
could be the Be star circumstellar disk, from which the neutron 
star accretes. 

The energy of the Fe fluorescence line does not vary. We 
obtained a weighted mean of 6.45±0.03 keV. The strength of the 
Fe gaussian line is strongly correlated with flux, pointing to an 
increasing fraction of reprocessed material by the circumstellar 
disk while the flux increases. The gaussian line width was 
almost constant during the outburst, and was eventually fixed 
at 0.5 keV in all the fits. The equivalent width (EW) of the line 
remained fairly constant during the outburst, with a weighted 
mean value of 0.12±0.04 keV, with the exception of the few 
observations in which A^h reached the highest values. In this 
case, also the iron line EW increased to higher values, around 
0.23 keV. 



3.2.1. Pulse-phase spectroscopy 

We performed spin phase-resolved spectroscopy for one of the 
longest pointings. The phase was divided into 10 bins, and for 
each bin a spectrum was extracted and fitted in the 3-20 keV 
energy range with the same components as the phase-averaged 
spectra. The cyclotron absorption feature was not included at 
this stage, since it lies outside the sensitivity range of the PCA. 
In the 3-20 keV range, its contribution is negligible. Five of the 
ten obtained spectra showed residuals at ~8 keV, consistent with 
an emission-like feature, which was not significant in the phase- 
averaged spectra. The physical origin of this feature is uncertain 
(ICoburn et al.ll2002t iRodes -Roca et al1l2009t IDoroshenko et al.l 
2010). This was fitted with an additional gaussian. Results are 
shown in Fig. [9] The modulation of the 3-20 keV flux reflects, 
as expected, the pulse profiles of Fig. [7] The spectral parameters 
clearly exhibit phase dependence: the photon index, the cutoff 
energy, A^h, and the strength of the iron line show larger values 
during the off-pulse phase, and lower values in correspondence 
to the on-pulse. The energy of the iron line did not show signifi- 
cant variation. 
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3.2.2. Cyclotron resonant scattering feature 

The cyclotron line was modeled using the GABS XSPEC com- 
ponent, a gaussian absorption line described by the following 
profile: 
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Fig. 10. Relation between the photon index and the QPO characteristic 
frequency (upper panel) and rms (lower panel). Different symbols mark 
the different phases of the outburst: rising (empty squares), peak (full 
circles), and decay (full triangles). 



where E c is the line energy in keV, cr is the line width in keV, 
and t is the optical depth at the line center. The obtained central 
energy of the line was ~60 keV, with some dispersion during the 
outburst. No harmonics were found in our energy range. This 
feature was significantly detected during the whole outburst, 
with the exception of the last four observations, corresponding 
to an X-ray flux lower than 2 x 10~ 9 erg citT 2 s _1 . 

In contrast to other X-ray pulsars , (see, for instance , 
lMiharaetal.ll 19981: iTsvgankov etal1l2006t IStaubert et aL"ll2007l) . 
there is no clear correlation between the central energy of 
the CRSF and the X-ray luminosity, although, on average, 
E c ~ 65 keV at low luminosities, and E c ~ 55 keV at high 
luminosities. Errors are quite large for this variation to be 
statistically significant. Nevertheless, the observed trend of 
higher energy at low luminosity agrees with the relation seen in 
other HMXB pulsars that show an anticorrelation between E c 
and luminosity. The value of ~60-65 keV is the highest energy 
measured for a cyclotron absorption feature of an X-ray pulsar. 
We thus confirm and extend to lower luminosities the findings 
by iDoroshenko et al.l d2010t) . who found a value of ~55 keV 
during the peak of the outburst. 

The presence of CRSFs provides a tool for direct measure- 
ment of the magnetic-field strength of accreting pulsars from the 
relation 



E c = 11.65i2X(l+z) -1 keV, 



(3) 



1(E) = exp(-(r/ V2^ cr) exp(-0.5((£ - E c )/cr) 2 )), 



(2) 



where B\2 is the magnetic field strength in units of 10 12 G 
and z is the gravitational re dshift in the line-forming region 
dWasserman & Shapirol 1983b . Assuming a gravitational redshift 
of z = 0.3 for a typical neutron-star mass of 1 .4 Mq and radius of 
10 km, we estimated the magnetic field to be 6.7xl0 12 G. Due to 
the high-energy CRSF, the value obtained for the magnetic field 
is among the highest measured for an X-ray pulsar. 

3.3. Correlations between spectral and timing characteristics 

Figures [10] and QT| show the relation between the photon index 
and the parameters of the QPO and the L\ component, respec- 
tively. The frequency of both components shows a negative cor- 
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relation with the photon index; the L\ rms shows a correlated 
behavior with T, while rmsQpo remains constant. The stronger 
correlation seen for L\ can be attributed to a wider range of vari- 
ability of F, which, in turn, results from a wider range of vari- 
ability in the X-ray flux. 



4. Discussion 

We have performed detailed spectral and timing analysis of 
the accreting pulsar 1 A 1 118-615 during an entire giant outburst. 



4. 1 . Aperiodic variability analysis 

We reported the discovery of a QPO in the power spectra of the 
system (Fig. [2}. This peaked component has a characteristic 
frequency of 0.08 Hz, and it is only observable in high-flux 
observations, down to Lx ~ 40% of the maximum observed 
luminosity. Including 1A 1118-615, QPOs have been detected 
in 17 accretion-powered high-magnetic field pulsars, com- 
prising many HM XBs and a fe w LMXBs, both transient and 
persistent (see lJames et alJuOlOl and references therein). QPOs 
in X-ray pulsars provide strong evidence for the presence of an 
accretion disk during giant outbursts because they are believed 
to be related to inhomogeneities in the inner disk dPaul & Raol 
1 1998b . While black-hole binaries and LMXBs show QPOs with 
frequencies from a few Hz to a few hundred Hz, high-magnetic 
field neutron stars only show low-frequency QPOs, in the range 
between 10 mHz and about 1 Hz. 

The presence of QPOs in the power spectral density of X-ray 
pulsars is in general explained with either the Magn etospheric 
Beat Frequency Model (BFM, lAfoar & Shahaml 1 19851) or th e 
Keplerian Frequency Model (KFM. Ivan der Klis etaD Il987l) . 
The main difference between the two models is the interpre- 
tation of the QPO. In the BFM, vqpo = Vk — v s , where Vk is 
the Keplerian frequency of the inner accretion disk, and v s is 
the spin frequency, while in the KFM, vqpo = vt. Given that 
in 1A 118-615 v s <K vqpo, it is very difficult to distinguish 
between the two models. A third mode l, the magnetic disk 
precession model (Shirakawa & Lai 2002) attributes mHz QPOs 
to warping/precession modes induced by magnetic torques near 



the inner edge of the accretion disk. None of the three models 
can explain the origin of QPOs in all the systems. 

If one assumes that QPOs are produced as a result of 
Keplerian motion of inhomogeneities in an accretion disk, 
then the longer implied timescales of the QPOs in HMXBs are 
somehow expected since the inner radius of the accretion disk 
must be larger than the magnetospheric radius Rm ~ 10 8 cm. 
The inner disk radius can be expressed as ro = (GM/47r 2 v^)^ 3 , 
where is the Keplerian rotation frequency at the inner edge. 
An 80 mHz QPO implies a disk radius of 8-9xl0 8 cm, that is, 
outside the magnetosphere, as models foresee. 

The QPO frequency can be employed to estimate the 
neutron-star magnetic field. The Alfven radius for a 1.4 Mq 
neutron star with radius of 10 km can be approximated to 
R M - 3 x 10 8 L3 ? 2/ Vto 7 cm dFranket alj|2002h . where L 37 is 
the X-ray luminosity in units of 10 37 erg s _1 , and //30 is the 
magnetic moment in units of 10 30 G cm 3 . Assuming that the 
QPO is formed at the inner accretion disk, that is, equating the 
inner accretion disk radius as obtained from the QPO frequency, 
ro, with Tm, we estimated a magnetic moment of 7.8 xlO 30 G 
cm 3 , which is equivalent to a magnetic field of 7.8 xlO 12 G. 

A positive correlation is expected between the QPO fre- 
quency center and the X-ray flux, because at larger mass 
accretion rate, corresponding to increasing flux, the size of the 
magnetosphere decreases, and the accretion disk is expected 
to extend closer to the neutron star. This implies shorter 
characteristic time-scales, which means higher frequencies. In 
the case of 1A 1118-615 during the outburst analyzed in this 
work, we have indeed detected a positive correlation between 
the QPO frequency and the 3-30 keV flux (Fig. 3). In contrast, 
the fractional rms amplitude does not show any correlation with 
flux. Note that the overall rms of the continuum remained fairly 
constant during the outburst. Thus, the disappearance of the 
QPO cannot be linked to any overall increase in power, but must 
be attributed to a decrease in the strength of the QPO itself. 

We found a negative correlation with energy above 4 keV 
(Fig. 4). No general trend has been identified for QPOs in 
HMXBs, with cases of correlation (as V0332+53, see lOu et al.l 
120051). anti-correlation (XTE J1858+034 KS 1947+300, see 
Muk heriee et al.l 120061: iJames et all [2010). or no dependence 
(EXO 2030+375. see lAngelini et aLlll989h . 

The characteristic frequencies of the two broad-band noise 
components increase as the flux increases (Fig. 5). This is 
particularly true for L\, the best constrained and always present 
component. The correlated behavior is at least detected up to 
flux ~6xl0~ 9 erg ctrT 2 s -1 (L x = 1.8 x 10 37 erg s" 1 , for a 
distance of 5 kpc). If one assumes that all the aperiodic phe- 
nomenology arises in the same region, the correlation observed 
between frequencies and X-ray flux is somehow expected, and 
can be explained in the same manner as for the QPO frequency, 
namely, due to the shrinking of the magnetosphere as flux 
increases. 

The fractional amplitude of variability of the L\ component 
decreases as the flux increases (Fig. 5). This trend was not seen 
in the QPO rms due to the limited flux range in which the QPO 
was detected. This is also true for the L2 component, although 
in this case other effects may be at work. Changes in the mass 
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accretion rate will affect first large radii, i.e., low frequencies. 
Large amplitude variations may then be suppressed as they 
propagate inward. For fluxes larger than ~ 7 x 1CT 9 erg crrT 2 
s _1 , the correlation of frequencies and variability amplitudes 
with flux is no more observed. We propose that this saturation 
corresponds to the minimum extent of the magnetosphere, or, 
equivalently, to the minimum inner accretion disk radius. 



4.2. Spectral analysis 

The most remarkable feature in the X-ray energy spectra of 
1A 1118-615 is a cyclotron line at ~60 keV, the highest value 
measured in a HMXB. This valu e is consistent with recent 
results by iDoroshenko et al.l d2010h . Before this discovery, the 
X-ray pulsar showing the hig hest CRSF energy w as A0535+26, 
with a value of 52 keV dPostnov etall 120081) . It has been 
pointed out that the energi es of cyclotron absorption features 
vary w ith X-ray luminosity (Mihara et al. 2004; Nakaiima et al.l 
2006). Such variations are interpreted in terms of the height 
of the line-forming region above the neutron star surface. 
Both positive and negative correlations between the cyclotron 
line energy and X-ray luminosity were observed in X-ray 
pulsars, with a prevalent anti-correlation in high-luminosity 
pulsars and correlation in low-luminosity ones. Prominent 
examples for such a variability are V0332+53 and 4U 01 15+63 
on one sid e, where the cyclotron line energy decreases with 
luminosity (Tsygan kov et al.l 120061; [Nakaiima et al. 2006), and 
Her X- 1 on the other side, where the line energy increases with 
luminosity (Staubert et aH 120071) . A counterexample is on the 
other hand represented by A0535+262, where the line energy 
does not show any kind of correlation with flux C Caballero et al.l 
[2007b . In the case of 1A 1118-615, no strict trend associated 
with flux was detected, although higher values for the central 
energy were measured at lower luminosity. The properties of 
cyclotron lines are directly related to the pulsar X-ray emission 
processes and the behavi or of matter in strong magnetic fields 
(Meszaros & Nagellll985h . Therefore, CRSFs carry a wealth of 
information about the environment in which they are formed. 
CRSFs provide the only direct measurement of the neutron star 
magnetic field, that we calculated to be ~6.7xl0 12 G. Given 
the order-of-magnitude calculation, this is compatible with 
the value of 7.8xl0 12 G, obtained independently from timing 
analysys. 

The photon index anti-correlates with flux (Fig. [8}, in 
contrast to what was found in other Be/XR Bs (like, for in- 
stance, EXO 2030+ 375, V0332+53, see IWilson et al.l [2008; 
IMowlavi et al.l 120061) , but similar to what was observed in 
SAX J2103.5+4545 dReig et al.l |2010|) . This trend implies a 
hardening of the spectrum as the flux increases. The hard emis- 
sion from X-ray pulsars is believed to originate in the accretion 
column, as a result of Comptonization of soft photons coming 
from the neutron star thermal mound by hig h-energy electrons 
in the accretion flow dBecker & W olff 2007). We speculate that 
as the accretion rate increases, the photons spend more time, 
on average, being upscattered in the flow before escaping from 
the accretion column, thus resulting in harder X -ray spectra 
dLanger & Rappaportlll982l;lBecker & Wolff1l2007l) . 

The measured hydrogen column density was almost constant 
(A^h at ~3xl0 22 crrT 2 ) with time, with the exception of a few 
observations at the beginning and at the and of the outburst, 
where it more than doubled its value (A^h at ~7xl0 22 cm -2 ). 



After obtaining very poor fits by constraining the value of A^h 
to 3xl0 22 crrT 2 for all pointings, we believe that the observed 
variation is a real phenomenon. We measured the duration of 
this 'Wh transition", from lower to higher values, obtaining 
~3 days in both cases, at the beginning and at the end of the 
outburst. The time interval with A^h confined at lower values 
is approximately 15 days, between the two transitions. This 
phenomenon is a remarkable spectral feature of 1A 1118-615 
and must arises from some change in the matter surrounding 
the neutron star during the outburst, possibly due to the Be 
companion activity. It is very unlikely that the observed 15 
days interval has an orbital ori gin since, acc ording to the pulse 
period/orbital period diagram dCorbetll 19861) . a 407 s P s would 
imply a 400-800 days P or t,. Also, the lack of Type-I and scarcity 
of Type-II outbursts seen in thi s system would contrast wit h 
such a short orbital period (see Okazaki & Negueruela 2001). 
More measurements by observatories with high sensitivity at 
low energies are needed to investigate this peculiar property. 
The hydrogen column density as obtained from our data was 
compared with the value converted from the visual extinction 
Ay. A difference of one order of magnitude points to a source of 
absorption that is local to the compact object, consistently with 
the on-going X-ray giant outburst. 

The center of the Fe Key-emission line was determined 
to be 6.45+0.03 keV during the entire outburst. This value 
is compatible with neutral or weakly ionized iron. The line 
equivalent width is approximately constant with flux, while 
the strength of the line grows linearly with flux (Fig. |S}. This 
result gives firm evidence that the iron line in fact arises from 
the fluorescence of continuum X-rays absorbed by cool matter 
surrounding the central source. High hydrogen column density 
and the presence of the Fe line during all the outburst attest to a 
large amount of absorbing material around the compact object. 
This is a further evidence of the extremely dense and extended 
disk around the Be companion, whi ch is the Be star d isplaying 
the strongest Balmer emission lines dMotch et al.lll988h . 



4.3. Pulse profile analysis 

We also studied the pulse profiles at different energies, finding a 
dependence of the pulse shape with the energy band (Fig. [7]). A 
double peak structure is observed below ~ 10 keV, while a single 
peak profile is detected above that energy. The pulsed fraction 
increases with energy from ~30% to ~65%, for energy above 
~10 keV, with a "pla t eau" a t lower energies, in agreement 
with Doroshenko et al. (2010), similar to what was found in 4U 
0115+63 dFerrigno et al.ll2009l) . Various types of pulse profiles 
are observed in X-ray pulsars. The profile energy dependence 
displayed by 1A 1118-615 is of the same kind as, for i nstance, 
4U 115+63 dTsygankov et all 120071) and Her X-l dNagasel 
1989). The presence of a second, weaker, peak is usually 
interpreted as the parti ally obscured emission of the column 
of the antipodal spot (Tsygankov et al. 2007). Pulse profile 
changes may have various origins, like changes in the spot 
shape and position, variations of the emission pattern, of the 
inner disk radius or the optical depth of the accretion stream. 

To investigate the dependence of the main spectral param- 
eters with the spin phase, we performed spin phase-resolved 
spectroscopy (Fig. [9j. Our results show a clear modulation of 
spectral parameters with the spin phase. In particular, spectra 
are softer at off-pulse phases, and harder at on-pulse phases. The 
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cutoff energy, the hydrogen column density and the intensity 
of the iron fluorescence line display the higher values at the 
off-pulse phase. We speculate that a hardening of the spectrum 
as the spin phase gets on-pulse may be due to the fact that the 
accretion column, and thus the Comptonized X-ray emission, 
are being observed from a more direct and somehow transparent 
line of sight. On the other hand, at the off-pulse phase, much 
more material interpose between the emission flow and the 
observer, resulting in softer spectra, higher hydrogen column 
density, and larger quantity of reprocessed iron, as the spectral 
parameters coherently indicate. 



4.4. Spectral/timing correlations 

Figure 1 1 shows the relationship between the photon index and 
the frequency and rms variability of the L\ noise component. 
This relationship results naturally from the correlations reported 
above of the photon index and the timing parameters with flux. 
However, Fig. 1 1 allows a quick comparison of the X-ray tem- 
poral variability as a function of spectral shape. This is the first 
time that such correlation is reported in an accretion-powered 
pulsar. Figure 10 shows the same relationship for the QPO. 
The correlation is much weaker than for L\, possibly due to the 
limited flux range where the QPO is observed. 

Figure 1 1 indicates that the source is more variable when it 
shows a soft spectrum, i.e., at low flux and, more importantly, 
that the physical region where the aperiodic variability origi- 
nates and the emission region are strongly linked to each other. 
The aperiodic variability is supposed to arise from a region 
outside the magnetosphere at R > 10 s cm, while the emitting 
region that originates the energy spectra lies in the accretion 
column, i.e., at distances close to the neutron star surface, 
namely 10 6 - 10 7 cm. Nevertheless, these two regions are 
somehow physically connected. The correlation between rms 
and photon index might be explained by simple Comptonization 
models, if one assumes that the source of variability is variations 
in the soft photon input. In this case, low-energy photons would 
retain most of the variability of the seed photon input since 
they did not spend much time in the Comptonization medium 
and have not undergone many scatterings. In contrast, the 
variability of high-energy photons would be smeared out since 
these photons have spent longer time in the accretion flow. The 
overall result would be higher variability in correspondence to a 
soft spectrum, as observed. 

This work demonstrates the need for a comprehensive ap- 
proach to the study of accretion-powered pulsars. The complex- 
ity of the models of accretion in these systems, mainly due to the 
strong coupling of the magnetic field and the accretion flow, cru- 
cially requires the employment of all kinds of techniques, multi- 
band photometry (colors), spectral and timing analysis. 
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